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Abstract 

O ■ 

The use of rapidity gaps is proposed as a measure of the spatial pattern of an 
r-H ■ event. When the event multiplicity is low, the gaps between neighboring particles 

carry far more information about an event than multiplicity spikes, which may occur 
very rarely. Two moments of the gap distrubiton are suggested for characterizing an 
■ event. The fluctuations of those moments from event to event are then quantified by an 

entropy-like measure, which serves to describe erraticity. We use ECOMB to simulate 
the exclusive rapidity distribution of each event, from which the erraticity measures 
£Sj ' are calculated. The dependences of those measures on the order of q of the moments 

provide single-parameter characterizations of erraticity. 

ON 
Oh! 



1 Introduction 



To study the properties of event-to-event fluctuations in multiparticle production, it is nec- 
essary to have an effective measure of the characteristics of the final state of an event. The 
totality of all the momenta of the produced particles constitutes a pattern. A useful measure 
of a pattern should not contain too much details, but enough to capture the essence that 
is likely to fluctuate from event-to-event. In previous papers we have used the normalized 
factorial moments, F q , as a measure for studying chaos in QCD jets Jl[, in classical nonlinear 
dynamics 0, and erraticity in soft production of particles [[|. We now consider two different 
moments in order to improve the analysis in problems where the use of F q is less effective. 

The horizontal factorial moments register the multiplicity fluctuation from bin to bin in 
an event. However, when the event multiplicity is low, and the bin size small, most bins have 
only one particle per bin, and the factorial moments fail to provide a good characterization of 
the event pattern. To overcome that deficiency, we shift our emphasis from bin multiplicities 
to rapidity gaps. It is intuitively obvious that the two quantities are complementary: the 
former counts how many particles fall into the same bin, while the latter measures how 
far apart neighboring particles are. Clearly, the former works better when there are many 
particles in an event, while the latter is more suitable when there are few particles. 

The search for a good measure of event-to-event fluctuations can be carried out only if 
we have an event generator that can be used for the exploration. To that end we shall use 
ECOMB [|J, which simulates soft production processes in hadronic collision. It is based on 
a reasonable modeling of the many-body dynamics in which the partons undergo successive 
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color mutation before hadronization. It is the only model capable of generating factorial mo- 
ments that agree with the intermittency data of NA22 ||. However, there is still freedom in 
the model for further adjustment. The aim of this paper is not to test ECOMB or to improve 
it. Despite its imperfections, it can nevertheless simulate events with sufficient dynamical 
fluctuations that deviate significantly from statistical fluctuations. That capability is what 
we utilize in our search for the desired measure. For that reason it is unnecessary for us to 
review here the dynamical content of ECOMB. After the experimental data are analyzed 
and the proposed measure determined, we can then return to the problem of modeling soft 
interaction when the new erraticity data will provide the guidance needed for an upgrading 
of ECOMB. 

2 The Problem 

Let us start by reviewing the factorial moments for multiparticle production ||. They are 
defined (for the gth order) by 

f q =(n(n-l)---(n-q + l)) , (1) 

where n is the multiplicity in a bin. Originally, the average in ([!]) is performed over all events 
for a fixed bin, which we now call vertical average. Later, the horizontal average is considered 
for a fixed event, where n in ([!]) is averaged over all bins. In either case if the probability 
distribution P n of n can be expressed as a convolution of the dynamical distribution D{y) 
and the statistical (Poissonian) distribution, i.e., 

P n = jdv^v n e- v D(v) , (2) 

then one obtains || 

f< = £^ p > = f*'" D M ■ (3) 

Since it is a simple moment of D(u), the statistical fluctuation is regarded as having been 
filtered out by f q . 

The above procedure of eliminating the statistical fluctuation fails either when the sum in 
(§) does not extend to infinity, or if that fluctuation cannot be represented by a Poissonian 
distribution as in (|2|). Both of these circumstances occur for horizontal analyses of low 
multiplicity events. There is nothing wrong with calculating f q according to (jH) for such 
events. The question is what one can use f q for. 

In the horizontal normalized factorial moments F q = f q /fl are used to characterize the 
spatial pattern of an event. Such a characterization clearly cannot convey all the details of an 
event; indeed, extensive details using many variables are not desirable for the quantification of 
event-to-event fluctuations. It is evident from (j^) that only bins with n > q can contribute 
to f q , but the positions of the contributing bins have no effect on f q . That deficiency is 
unimportant when many bins contribute. However, when the event multiplicity N is low 
and the number of bins M is high, so that the average bin multiplicity n = N/M <C 1, then 
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it is only by large fluctuations that a bin may have n > q, whether they are dynamical or 
statistical in nature. Since f q is insensitive to where the few contributing bins are located, 
there is very little information about an event that is registered in F q . In |7] it is shown that 
in the framework of a simple model F q are dominated by statistical fluctuations when N is 
small, but they reveal the dynamical fluctuations when N is large. 

The aim of this paper is to find an alternative to F q that can effectively characterize the 
spatial pattern of an event, even when the event multiplicity is low. 

3 The Solution 

From ([!]) we see that f q receives a contribution from a bin in which n > q, but ignores where 
it is located. In other words f q is sensitive to the local height of the rapidity distribution in 
an event, not to the spatial arrangement in rapidity. When N is low and M is high, many 
bins are empty. To have a bin with n > q means that even more bins than average would 
have to be empty. It then seems clear that the complementary information accompanying 
rapidity spikes is the rapidity gaps. When N is high, rapidity gaps are generally not very 
informative; however, when N is low, they characterize an event better than counting spikes. 
In the following we shall develop two methods based on measuring the rapidity gaps. 

Since particle momenta can be measured accurately, there is no need to consider discrete 
bins in the rapidity space. Thus we shall work in the continuum. Moreover, the advantage of 
working in the cumulative variable X has long been recognized || ||, and we shall continue 
to use the X variable, as in M ( though not explicitly stated there in the first reference). 
The definition of X is 

X(y)= [" p{y')dy'/ phftdtf , (4) 

where p {%)') is the single-particle inclusive rapidity distribution and 2/ m in(max) is the minimum 
(maximum) value of y. Thus the accessible range of y is mapped to X between and 1, and 
the density of particles in X, dn/dX, is uniform. 

Consider an event with N particles, labeled by i — 1, • • • , N, located in the X space at 
Xj, ordered from the left to the right. Let us now define the distance between neighboring 
particles by 

Xi = X i+1 -Xi, i = 0,---,N , (5) 

with Xq = and X i+1 = 1 being the boundaries of the X space. Every event e is thus 
characterized by a set S e of N + 1 number: S e = {xi\i = 0, • • • , N}, which clearly satisfy 

f> = l . (6) 

i=0 

We refer to these numbers loosely as "rapidity" gaps. 

For any given event S e contains more information than P n , which is the bin-multiplicity 
distribution for that event; in fact, P n can be determined from S e , but not in reverse. To 
study the fluctuation of S e from event-to-event is the most that one can do; indeed, too 
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much information is conveyed by S e . For economy and efficiency in codifying the informa- 
tion we consider moments of Xi that emphasize large rapidity gaps. As mentioned earlier, 
concomitant to the clustering of particles that results in spikes in the rapidity distribution is 
the existence of large gaps. Thus moments that emphasize large Xi convey similar informa- 
tion about an event as do the moments that emphasize the high-n tail of P n . However, the 
factorial moments of P n suffer the defects discussed in Sec. 2 that are absent in the moments 
of Xt. The issue of statistical fluctuations has to be addressed separately. 
Let us then define for each event 



G 



1 A 



N + . „ 



r X>? > (7) 



Despite the similarity in notation, these moments bear no relationship to the G-moments 
considered earlier |1(J. It is clear from (|6|) and fl7|) that 



Go = 1 and G\ 



N + 1 



At higher q, G q are progressively smaller, but are increasingly more dominated by the large 
Xi components in S e . A set of G q for q ranging up to 5 or 6 is sufficient to characterize 
an event, better than S e itself in the sense that G q can be compared from event-to-event, 
whereas S e cannot be so compared due to the fluctuations in N. 
If we define the gap distribution by 

1 N 

9{x) = jt— tJ2 S ( x ~ X i) > ( 9 ) 

then the G moments are 

G q = ( 1 dxx q g{x) . (10) 





This form may become more convenient in some situations. 

Since G q fluctuates from event to event, we can determine a distribution P(G q ) of G q 
after many events. It is the shape of P{G q ) that characterizes the nature of the event-to- 
event fluctuations of the gap distribution, and therefore of the spatial pattern of an event. 
Again, we can describe P(G q ) by its moments 

c m = JfTX Ge q ) P = J dG q G p q P(G q ) , (11) 

where e labels an event and M is the total number of events. Since we need not consider 
bins in x, G q is a number for each event without statistical error. Thus calculating the 
pth moment does not compound statistical errors. Although one can consider a range of p 
moments, we shall focus only on the derivative at p — 1 in the following. 

Since C\ q = (G q ) is the mean that gives no information on the degree of fluctuation, the 
derivative at p — 1 convey the broadest information on P(G q ). We have 



s q = -j- p C P , q \ p=1 = -{G q lnG q ) , (12) 
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where (• • •) stands for average over all events. The quantities s q are our new measures of 
erraticity in terms of rapidity gaps. Since G q is not a probability distribution, s q is not an 
entropy function, despite its appearance. 

Unlike the factorial moments, G q does not filter out statistical fluctuations. At low 
multiplicities F q fails to be effective in that filtering anyway, as discussed in Sec. 2, so it is 
at no great loss to consider G q . However, we can have an estimate of how much s q stands 
out above the statistical fluctuation by first calculating 

sf = -(GflnGf) , (13) 

where G s q l is determined from (|ToD by using only the statistical distribution of the gaps, 
g st (x), i.e., when all N particles in an event are distributed randomly in X space. Then we 
take the ratio 

S q = s q /sf , (14) 

and examine how much S q deviates from 1. S q will be the first erraticity measure that we 
shall calculate in the next section. 

Since our interest is in the deviation of G q from (G q ), a measure of that deviation is 

^-(<t>'"t>) ■ (15) 

which clearly would be zero if G q never deviates from {G q ). We can further normalize s s q l by 
the statistical-only contribution s s q l and define 

S q = s q /~sf . (16) 

Whether this is a better quantity to represent erraticity will be examined quantitatively in 
the next section. 

The moments G q are not the only ones that can characterize the rapidity-gap distribution. 
In fact, since Xi < 1, G q are usually <C 1, and the statistical errors on S q and S q turn out to 
be quite large, though not so large as to render the measures ineffective. We now consider 
a different type of moments that also emphasize the large gaps. Define for an event with N 
particles 



H, 



1 v 



s 



where Xi is as given in (f|). These moments also receive dominant contribution from large 
do G q , but H q can become ^> 1. In terms of g{x) we have 

H q = I dx(l-x)- q g(x) , (18) 
Jo 

where g{x) must vanish sufficiently fast as x — > 1 to safeguard the integrability of 
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We can substitute H q for G q in all of the foregoing considerations. In particular, we can 
define 



a q = {H q lnH q ) , (19) 

' (20) 

S, = ^ and E 9 = 4 (21) 
of af 



as new measures of erraticity. The only nontrivial point to remark on concerns the event 
average. 

For each event H q depends implicitly on the event multiplicity N. If V n is the multiplicity 
distribution, then the average (H q ) is given by 

oo 

(Hq) = E H g (N)V N , (22) 

N=q+l 

where H q (N) is the mean H q after averaging over all events with N particles in each event. 
Note that the sum in fl22|) begins at N = q + 1, not 0. To see this subtle point, let us start 
with the statistical average for which we can make precise calculations. In the Appendix 
we show that the probability distribution Pn(x) of the gap distance x, after sampling with 
sufficiently many events, each with N randomly distributed particles in the X space, is 

p§(x) =N(l-x) N ~ 1 . (23) 

Thus it follows that 

r 1 N 

H s q \N)= dx(l-x)-«p s *(x) = - . (24) 

Jo N — q 

Evidently, iV must be greater than q to ensure convergence. If for statistical calculation we 
require n > q + 1, then we make the same requirement for the general problem in (^), so 
that <j q and a s q l in ( pip are calculated on the same basis. 

4 Results 

We have applied ECOMB j4|, upgraded by [[□]], to calculate the rapidity distribution for each 
event. From that we compute the gap distribution g(x) in the X space. After simulating 10 6 
events at y/s = 20 GeV, our result for S q is shown in Fig. 1. The error bars are determined 
by using the conventional method. The straight line in Fig. 1 is a linear fit of the central 
points. Evidently, the result indicates a power-law behavior in q for q > 2 

S q oc q a , a = 0.156 . (25) 
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The fact that S q deviates unambiguously from 1 implies that it is a statistically significant 
measure of erraticity in multiparticle production. At y/s = 20 GeV the average charge 
multiplicity is only 8.5, which is low enough to cause problems for the factorial moments F q , 
but our use of the gap moments G q evidently encounters no similar difficulty 

We next consider S q defined in ([15]) and ([RJ) . The result is that S q is nearly independent 
of q, as shown in Fig. 2. More precisely, we obtain 

S q = 0.96 ± 0.03 . (26) 

We regard this result as indicative of the inadequacy of S q as a measure of erraticity, since 
S q is almost consistent with 1. 

Turning to the H q moments, we show in Fig. 3 in semilog plot the dependence of S g on 
q. Evidently, a very good linear fit is obtained, yielding 

S 9 oc e Pq , 0= 0.28 . (27) 

In the same figure we also show S 9 . Although the error bars are larger, an exponential 
behavior 

t q oc e h , (3 = 0.25 , (28) 

can nevertheless be identified. Note that E ? is much farther from 1 than S q . Since Eg has 
less statistical error than E 9 , it is more preferred. Hereafter we shall discard S q and E^ from 
any further consideration. 

We now examine the dependence on cm. energy. The higher multiplicities at higher s 
will decrease the average gap (x) and the corresponding moments G q and H q . We calculate 
the effects on S q and E g at y/s = 200 GeV. The results are shown in Figs. 4 and 5, where 



the values for y/s = 20 GeV are reproduced for comparison. The power law (|25|) and the 



exponential behavior (|27|) persist; the corresponding parameters are 

a = 0.133 , /? = 0.108 at ^ = 200 GeV . (29) 

Whereas a has changed little, j3 has decreased significantly. The variability of (3 makes 
it a more sensitive measure of erraticity, although the stability of a may nevertheless be 
interesting and useful. Only the analysis of the experimental data will reveal which one 
between S q and E g is better in quantifying erraticity. It can also turn out that both are 
good. 



5 Conclusion 

We have proposed the moments G q and H q as measures of spatial patterns in terms of rapidity 
gaps. We then showed that the entropy-like quantities S q and E 9 deviate sufficiently from 1 
with small enough statistical errors to serve as effective measures of erraticity, i.e., event-to- 
event fluctuations. In the framework of an soft hadronic interaction event generator ECOMB 
we have obtained the behaviors S q oc q a and E 9 oc e^ q . The precise forms of these results are 
unimportant from the point of view of the search for an experimental measure to quantify 
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erraticity. We offer both S q and S g as our findings. On the other hand, from the point of 
view of using erraticity to test event generators, then the forms of our results for S q and S g 
are pertinent, and the values a = 0.156 and (3 = 0.28 are useful for comparison with the soft 
production data. Analysis of the data, especially those of NA22, to determine S q and S q is 
therefore urged. The experimental values of a and j3 will either eliminate wrong models or 
provide crucial guidance to the improvement of the correct models. 

The extension of this approach to other collision processes is obviously the next step. For 
heavy-ion collisions the multiplicities will be too high for any interesting study in rapidity 
gaps, unless one focuses on more rarely produced particles, such as J/ip. At RHIC when only 
pp collisions are studied, our result for y/s = 200 GeV can be tested. For nuclear collisions 
very narrow Apx selection must be made to render the rapidity gap analysis meaningful. 

A natural direction of generalization is, of course, to higher dimensional analysis. One- 
dimensional gaps should be generalized to two-dimensional voids, which is more difficult to 
define if the use of bins is to be avoided. When a good measure is found, not only can it 
be employed as an alternative to the multiplicity analysis in the lego plot, useful application 
can no doubt be found also in the study of galactic clustering in astrophysics. Finally, in 
view of the abundance of experimental data and the variety of event generators for e + e~ 
annihilation, a generalization to the multidimensional variables suitable for such problems 
will be a fruitful direction to pursue. 
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Appendix 

We derive in this Appendix the probability distribution of gaps in the purely statistical 
case. 

The gap distribution g(x) defined in @ and @ can be more elaborately written as 
g e (x; Xi, ■ ■ ■ , Xn) for the eth event with iV particles located at X%, ■ ■ ■ ,Xn- It has the 
normalization 

/ dxg e (x;Xi, • ■ ■ ,X N ) = 1 . (30) 
J o 

If the value of Xj is randomly selected in the interval < Xj < 1, then after a large number 
of events J\T with N particles in each, probability distribution in x is 

1 N 

Pn( x ) = -jr?z29e{x]X 1 ,---,X N ) 

J * e=l 

= N\ f 1 dX 1 f 1 dX 2 ■ ■ ■ f 1 dX N g e (x; X ir --,X N ) , (31) 

in which the primitive distributions of the individual Xj values that should appear inside 
the integral have been set equal to 1 due to the statistical nature of their occurrences. 

Consider a specific gap y = Xj — Xi, where j — i + 1. Then for all the multiple integrals 
at and before Xi, we may reverse the order of integration and obtain 

dX l C dX 2 --- f 1 dXi = f 1 dXi \ Xl dX^ ■■■ [ X2 dX 1 = f 1 dXi—^—Xt 1 . (32) 
o Jx 1 Jx^ Jo Jo Jo Jo {i — ly. 

For all the multiple integrals at and after Xj, we have 

^■■■JL^^why^-^' 1 ■ (33) 

Substituting these and (|5|) and (0) into (0), we get 

/VI N r l r*i X i-1 (\ — Y 

= xTTgio^io dx *TT^)\ (N-j)i 5{x ~ v) ■ (34) 



i=0 • 

\N-1 



The integral over Xj yields (1 — y) 1 /{N — 1)!. Thus the final result is 

p^(x) = X(l-x) 7V " 1 . (35) 
This behavior is verified by numerical simulation. 



9 



References 

[1] Z. Cao and R. C. Hwa, Phys. Rev. Lett. 75, 1268 (1995); Phys. Rev. D 53, 6608 
(1996); ibid 54, 6674 (1996). 

[2] Z. Cao and R. C. Hwa, Phys. Rev. E 56, 326 (1997). 

[3] Z. Cao and R. C. Hwa, |hep-ph/990125l| Phys. Rev. D (to be published). 

[4] Z. Cao and R. C. Hwa, Phys. Rev. D 59, 114023 (1999). 

[5] I. V. Ajineko et al. (NA22), Phys. Lett. B 222, 306 (1989); 235, 373 (1990). 

[6] A. Bialas and R. Peschanski, Nucl. Phys. B 273, 703 (1986); 308, 857 (1988). 

[7] J. Fu, Y. Wu, and L. Liu, |hep-ph/99032T7| 

[8] A. Bialas and M. Gardzicki, Phys. Lett. B 252, 483 (1990) 

[9] E. A. DeWolf, I. M. Dremin, and W. Kittel, Phys. Rep. 270, 1 (1996). 

[10] R. C. Hwa, Phys. Rev. D 41, 1456 (1990); I. Derado, R. C. Hwa, G. Jancso, and N. 
Schmitz, Phys. Lett. B 283, 151 (1992). 

[11] R. C. Hwa and Y. Wu, Phys. Rev. D 60, 097501 (1999). 



Figure Captions 

Fig. 1 lnS q vs q as determined by ECOMB. The solid line is the best fit of the central 
points. 

Fig. 2 S g vs q with the same comments as in Fig. 1. 

Fig. 3 CnL q and lnL q vs q with the same comments as in Fig. 1. 

Fig. 4 inS q vs q at two different energies. 

Fig. 5 tnL q vs q at two different energies. 
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